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1. INTRODUCTION 


Image interpretation is an important method for acquiring train- 
ing data to classify Landsat Images in the Large Area Crop 
Inventory Experiment (LACIE) . Interpreting a scene for classi- 
fication requires that training fields containing statistically 
representative samples of all spectral signatures in the given 
scene be selected and correctly labeled. When multiple passes 
over a scene are to be interpreted, the variation of the spectral 
signatures in a multitemporal sense makes it difficult to select 
and identify all of the variety of signatures in the scene. To 
address this problem in multitemporal imaga interpretation of 
LACIE segments, the principal component (PCOMP) cluster images 
are introduced as an image interpretation aid. 

Some analyses of PCOMP transformed Landsat data have been per- 
formed (refs. 1,2). In this report, the PCOMP transformation 
was applied to some multitemporal clustered LACIE data. The 
object was to find whether PCOMP cluster images were helpful in 
selecting all multitemporal signatures and how much information 
(in terms of probability of misclas:' Ification) was lost when 
the three most significant PCOMP channels were used in classi- 
fying the PCOMP cluster image. 

A cluster image was generated by clustering the picture elements 
(pixels) in the segment and then replacing each pixel by the 
mean of the cluster to which it was assigned. Cluster keys were 
also generated as part of the cluster image (ref. 3). The keys 
were ordered according to the Kauth greenness number (ref. 4). 

A PCOMP cluster image was generated by applying the PCOMP trans- 
formation to the cluster image and the cluster keys. The PCOMP 
transformation was :he matrix of eigenvectors obtained from the 

mixture covariance matrix 1^. A cluster covariance matrix 

m c 

also was calculated using the cluster means in the cluster image. 



The purpose was to have an alternative method of obtaining the 
PCOMP transformation. It was found that E was not significantly 

c 

different from E ; hence, all PCOMP transformations in this 

in 

report were calculated from E„, 

m 

If m acquisitions from a LACIE segment were clustered in n chan- 
nels, where n « 4m, an n-dimensional cluster image (with 
n-dimensional cluster keys) would be generated. Notts that m 
color infrared (CIR) film products of the cluster image could be 
generated using a production film converter (PFC) . The first 
film product would use channels 1, 2, and 4. The second film 
product would use channels 5, 6, and 8; whereas the mth film 
product would use channels (n - 4) , (n - 3) , and n. A color film 
product of the segment PCOMP cluster image would be generated by 
the PFC using the first three (which are the most significant) 
channels of the image. CIR film products of the segment cluster 
image and PCOMP cluster image would also show color cluster keys. 

2. TEST AND EVALUATION PROCEDURES 

2.1 TEST PROCEDURE 

The following procedural steps were carried out in performing 
the test: 

a. For every image, the PCOMP transformation was calculated 
using E„, where E was the mixture covariance matrix cal- 
culated from all pixels in the image. 

b. For every image, a cluster image was generated using stand- 
ard clustering parameters. 

c. Using the PCOMP transformation, a PCOMP cluster image was 
generated . 

d. CIR film products were generated from the cluster image. 

e. One color film product was generated from the PCOMP cluster 
image . 
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f. CIR film products of the image were also generated. 

g. Using products from steps e, t, and g, along with the image 
ancillary data, color cluster keys were labeled as wheat or 
nonwheat. 

h. The initial wheat proportion estimate and nonwheat propor- 
tion estimate ^ computed. 

i. Using the labeled color keys in the PCOMP cluster image film 
product, representative wheat and nonwheat training fields 
were selected in proportion to and q^^. 

j. Ground truth (GT) test fields were selected for each segment 
using GT maps. 

k. An Earth Resources Interactive Processing System (ERIPS) 
batch job was submitted to classify each segment twice using 
the original features. One classification was performed 
using the training i'ields selected in step i with equal 
wheat and nonwheat a priori probabilities. This classifi- 
cation will be referred to as the PCOMP f ielde/equal a priorio. 
The second classification was performed using the training 
fields selected in step i with q^^ and q^^ obtained from step h, 
as the wheat and nonwheat a priori probabilities, respectively. 
This classification will be referred to as PCOMP fielda/uaer- 
input a prioria, 

£. The ERIPS was used to classify each PCOMP cluster image 

interactively using the first three PCOMP features. This 

(3) 

classification will be referred to as PCOMP fielda/equal 
a prioria. 

A flow chart of the test procedure is shown in figure 1. 
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Figure 1.— Flow chart for test procedure. 


4 















2.2 EVALUATION PROCEDURE 


a. Classification with original features: 

(1) The average probability of correct classification (PCC.„„) 

cl Vy 

over all test segments was calculated using PCOMP fields/ 
equal a prior ie procedures and compared with the PCC^^,q 
values obtained using the LACIE-selected training fields 
standard procedures. 

(2) The mean square error (MSE) in wheat proportions was 
calculated using PCOMP fielde/equal a priorie procedures 
and compared to the MSE values obtained using the LACIE- 
selected training fields standard procedures. The GT 
wheat proportions were used as a basis for calculating 
the MSE. 

(3) The average training field accuracy over all test seg- 
ments was calculated using PCOMP fielde/equal a priorie 
procedures and using the LACIE-selected training fields 
standard procedures. 

b. Classification with the first three PCOMP features: 

( 3 ) 

(1) The PCC^^g was calculated using PCOMP fielde/equal 

a priorie procedures and compared with the LACIE PCC^^^ 
and the PCOMP fielde/equal a priorie PCC^^^ values. 

(2) The MSE in wheat proportions was calculated using 
PCOMP^^^ fielde/equal a priorie procedures and compared 
with the LACIE and the PCOMP fielde/equal a priorie 
MSB's. 


3. DESCRIPTION OF THE DATA USED 

Most of the LACIE sample segments used in this test contained 
GT information. A description of those sample segments is shown 
in table I. Sample segments 1270, 1087, 1576, 1046, 1224, and 
1882 are blind sites. Sample segments 1976 and 1988 are intensive 
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TABLE I.- SAMPLE SEGMENT DESCRIPTION 


Acquisitions processed Location 


Segment 


(a) 

(county/state) 

1270 

76072, 

76162 

Throckmorton/Texas 

1087 

76001, 

76172 

Crosby /Texas 

1576 

76162, 

76198 

Lancaster/Nebraska 

1046 

76127, 

76190 

Beaver/Oklahoma 

1149 

76099, 

76189 

Harrison/Indiana 

1224 

76080, 

76162 

Dewey/Oklahoma 

1882 

76090, 

76153, 76189 

Lincoln/Kansas 

1976 

76136, 

76208, 76226 

Franklin/Idaho 

1988 

76038, 

76109, 76127, 76164 

Finney/Kansas 


^Numbers represent Julian dates; e.g, , 76072 is the 72nd 
day of 1976. 
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test sites (ITS's). Sample segment 1149 Is an operational LACIE 
segment for which no GT information is available. 

Each acquisition shown in table T represents 4 channels; hence, 

6 segments were processed in 8 channels, 2 segments were proc- 
essed in 12 channels, and 1 segment (1968) was processed in 
16 channels. 


4 . RESULTS 

The results of the test which is the subject of this report are 
presented in tabular form. Tables include the classification 
performance for each segment using the various procedures. The 
classification performance is specified in terms of training 
field accuracy; test field accuracy, which is the probability 
of correct classification (PCC) ; and wheat proportions. When 
equal a priori probabilities were used in classification, the 
training- and test-field accuracies were computed as the arith- 
metic average of the wheatfield and the nonwheatf ield accuracies. 
When LACIE procedures were used in processing the data set, 
equal a prioria were used in classification. When user-input 
a priori probabilities q^ and q^^ were used, the field accuracy A 
was computed as: 

* ’ * Vn 

and Aj^ was the nonwheat- 
field accuracy. It is worth mentioning at this point that the 
LACIE test-field accuracies for all blind sites were computed 
based on analyst-selected test fields, whereas the test-field 
accuracies using other procedures were based on fields selected 
from the GT maps. Analyst-selected test fields tended to show 
greater accuracy. 


where A^ was the wheatfield accuracy 
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The results of the test conducted are presented Ln tables II, 

III, and IV. Table II presents the training- and test-field 
accuracies for the data set using LACIE, PCOMP fielde /equal 
a priorie, and PCOMP^^^ fielde/equal a priorie procedures. It 
is noteworthy that (1) segment 1087 was processed manually in 
LACIE (i.e., the number of wheat pixels in the segment were 
hand counted) and (2) segments 1149 and 1988 were not processed 
multi temporally in LACIE. Table II shows that the LACIE and 
PCOMP fields /equal a priorie procedures yielded close training- 
and test-field accuracies. Table II also shows that the 
PCOMP fielde /equal a priorie procedure yields remarkably 

close accuracies tc those obtained when the PCOMP fielde /equal 
a priorie procedure was used. This shows that, in terms of 
PCC, no significant loss of information was noted when the 
first three PCOMP channels were used in classifying the data 
set used in this test. 

Table III presents the wheat proportions in the data set using 
GT and LACIE, PCOMP f ielde /equal a priorie, and PCOMP^^^ fields/ 
equal a priorie procedures. Table IV shows the effects of 
using the estimated a priori probabilities on classification 
performance. 

A sunimary of tables II, III, and IV is given in table V. For 
all the procedures used in this test, table V shows the average 
training-field accuracy, the PCC_„_, and the MSE in wheat 
proportions. The procedures displayed higher MSE in wheat pro- 
portions than the LACIE MSE, primarily because the wheat propor- 
tions in segments 1087 and 1576 were highly overestimated using 
the three PCOMP procedures. Table V indicates that the PCOMP 
fielde/ueer-input a priorie procedure was useful in the multi- 
temporal processing of LACIE segments. As shown also in table V, 
the PCOMP^ fields/ equal a priorie procedure displayed a good 
average classification performance in comparison with the rest 
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TABLE II.- TRAINING- AND TEST-FIELD ACCURACIES FOR THE DATA SET USING LACIE, PCOMP FIELDS/ 
EQUAL A PRTORIS, AND PCOMP^ FIELDS /EQUAL A PRIORIS PROCEDURES 
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TABLE III.- WHEAT PROPORTIONS IN THE DATA SET USING GROUND TRUTH AND LACIE 
PCOMP FIELDS/EQUAL A PFIOFISt AND PCOMP FIELDS/ 

EQUAL A PFIOFIS PROCEDURES 
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This segment was manually processed by the LACIE. 

’Processed by the LACIE using the Julian acquisition date shown in 





TABLE V.— SUMMARY OF RESULTS 
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of the procedures. This implies that most of the information 
contained in the multitemporal images for this test can be con- 
centrated in the first three PCOMP's, This test indicates that 
the PCOMP cluster images should be useful in the multitemporal 
processing of LACIE segments. 

5. CONCLUSIONS AND RECOMMENDATIONS 


Based on the results obtained in this test, the following conclu- 
sions were reached: 

a. Color film products of the first three PCOMP channels in the 
PCOMP cluster images are helpful interpretive aids for multi- 
temporal processing of LACIE segments. 

b. A procedure such as PCOMP f ielde/equal a priorie yielded 
comparable results to LACIE-selected training fields standard 
procedures. The PCOMP fields /user-input a prioris procedure 
displayed slightly better results than the PCOMP fields/equal 
a prioris procedure. 

c. It was found that most of the information (in terms of PCC) 
was preserved when the first three PCOMP channels were used 
in classifying the multipass PCOMP cluster images. 

It is recommended that the LACIE analysts be provided and become 
familiar with a color film of the PCOMP cluster image when multi- 
temporal processing of the LACIE image is required. 
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